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Abstract 

Background: MicroRNAs (miRNAs) play a fundamental role in the regulation of gene expression by translational 
repression or target mRNA degradation. Regulatory elements in miRNA promoters are less well studied, but may 
reveal a link between their expression and a specific cell type. 

Results: To explore this link in myeloid cells, mlRNA expression profiles were generated from monocytes and 
dendritic cells (DCs). Differences in mlRNA expression among monocytes, DCs and their stimulated progeny were 
observed. Furthermore, putative promoter regions of miRNAs that are significantly up-regulated in DCs were 
screened for Transcription Factor Binding Sites (TFBSs) based on TFBS motif matching score, the degree to which 
those TFBSs are over-represented in the promoters of the up-regulated miRNAs, and the extent of conservation of 
the TFBSs in mammals. 

Conclusions: Analysis of evolutionarily conserved TFBSs in DC promoters revealed preferential clustering of sites 
within 500 bp upstream of the precursor miRNAs and that many mRNAs of cognate TFs of the conserved TFBSs 
were indeed expressed in the DCs. Taken together, our data provide evidence that selected miRNAs expressed in 
DCs have evolutionarily conserved TFBSs relevant to DC biology in their promoters. 



Background 

In recent years, microRNAs (miRNAs) have taken center 
stage, as they are key regulators of gene expression at the 
post-transcriptional level, and play a fundamental role in 
a wide variety of biological processes, such as cell growth, 
development and several pathological conditions [1-3]. 
MicroRNAs are small, ~22 nt long, single-stranded mole- 
cules which, when complexed with an RNA-induced 
silencing complex (RISC), are able to form a complemen- 
tary double-stranded RNA structure by hybridizing to 
the 3' untranslated region of target transcripts, and inhi- 
bit translation of their cognate mRNA and/or promote 
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their degradation [4]. MicroRNAs have an established 
role in hematopoietic development and immunity. For 
example, forced expression of miR-181 in hematopoietic 
progenitors leads to an increase in the number of B cells 
[5], whereas it sets T cell receptor signaling thresholds by 
targeting negative regulators [6]. MicroRNA- 146a is up- 
regulated during toll-like receptor (TLR) signaling and 
targets TNF receptor-associated factor 6 (TRAF6) and 
IL-1 receptor-associated kinase 1 (IRAKI) [7], thereby 
serving in a negative feedback loop. Moreover, miR-155 
is also up-regulated during TLR and TNF signaling [8], 
and is required for normal immune function [9-14]. 

Although great strides have been made towards under- 
standing the biogenesis of miRNAs [4] and the identifica- 
tion of mRNA targets [15], their own expression is one of 
the least understood aspects. They are transcribed by 
RNA polymerase 11 [16] or RNA polymerase III [17]. In 
addition, approximately 80% of miRNAs are located in 
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introns of protein coding genes, but at least one third is 
believed to be transcribed independently from their host 
gene [4,18-20], whereas recent data suggest that most, if 
not all, intronic miRNAs contain putative promoters 
independent of their host gene [21]. In fact, it is now 
believed that once physically accessible, a gene is regu- 
lated by transcription factors that bind to their cognate 
transcription factor binding site (TFBS) in its promoter. 
Usually, there is more than one TFBS per gene, allowing 
combinations of transcription factors to elicit gene tran- 
scription. This phenomenon has been predicted for 
instance in Plasmodium falciparum, a parasite with a 
dearth of transcription-associated factors [22-24] and has 
been experimentally validated in other eukaryotic promo- 
ters [25,26]. 

Myeloid dendritic cells (DCs) and monocytes arise from 
a common monocyte/dendritic cell progenitor [27]. In 
vitro, DCs can be generated from blood-derived mono- 
cytes when cultured in the presence of the cytokines inter- 
leukin 4 (IL4) and granulocyte/macrophage colony 
stimulating factor (GM-CSF) [28]. DCs play an important 
role in innate immunity and the initiation of adaptive 
immune responses. They capture foreign antigens in per- 
ipheral tissues, migrate to the T-cell areas of secondary 
lymphoid organs and present these antigens to T- and B- 
cells. Depending on the extracellular signals they receive, 
they either induce tolerance in the steady state (tolero- 
genic DCs), or an inflammatory response in the presence 
of pathogen-associated patterns (PAMPs) or inflammatory 
cytokines (activated or mature DCs) [29,30]. As a conse- 
quence, DCs have gained considerable interest as vaccine 
adjuvants and are currently exploited in the treatment of 
cancer after loading with tumor-cell derived antigens 
[31,32]. 

In order to gain insight in miRNAs that may regulate 
DC development and behaviour, expression profiles of 157 
miRNAs were obtained from monocytes and DCs under 
inflammatory and tolerizing conditions. We show that 
DCs express a wide variety of miRNAs, some of which are 
differentially regulated during DC development and 
maturation. We predicted several target genes for these 
miRNAs, as well as binding sites for transcription factors 
in the putative promoter regions of these miRNAs. 
Furthermore, we show that by also taking evolutionary 
conservation [33,34] of the identified TFBSs into account, 
binding sites were found to preferentially cluster within 
500 bp upstream of the pre-miRNAs. Also, the fraction of 
conserved TFBSs for which the cognate transcription fac- 
tors are expressed in DCs increases with the number of 
miRNA promoters that contain these TFBSs. Taken 
together, the data described here provide evidence that the 
promoter regions of the miRNAs expressed in myeloid 
DCs contain binding sites for motifs of transcription fac- 
tors that are relevant to myeloid cell biology. This may 



help expand the understanding of the molecular mechan- 
isms underlying DC biology and development. 

Methods 

Isolation, culture and characterization of monocytes and 
DCs 

Peripheral blood mononuclear cells (PBMCs) were iso- 
lated as previously described [28]. Monocytes were iso- 
lated from PBMCs using CD14+ selection and the 
AutoMACS technology (Miltenyi Biotec, Bergisch Glad- 
bach, Germany) following the manufacturer's directions. 
Human monocyte-derived DCs were generated using 
GM-CSF and IL-4 and matured as described previously. 
Purity and maturation of monocytes and DCs were 
assessed by means of FACS analysis as described pre- 
viously [28]. Mixed lymphocyte reactions were essentially 
carried out as described elsewhere [35]. Briefly, they were 
co-cultured with peripheral blood lymphocytes in various 
dilutions for 2 to 4 days. Then tritiated thymidine (1 |iCi/ 
well; MP Biomedicals) was added to the cell cultures and 
incorporation was measured after 16 hr. Enzyme-linked 
immunosorbent assays (ELISA) were performed to assess 
the secretion of the inflammatory cytokines TNFa and 
IL8, using BD OptEIA (BD Biosciences, San Jose, CA) kits 
following the manufacturer's recommendations. 

Isolation of total RNA and microRNA-specific reverse 
transcription and quantitative PCR 

Total RNA was isolated using TRIzol (Invitrogen, Carls- 
bad, CA), following the manufacturer's recommendation. 
Quantity and purity were determined spectrophotometri- 
cally. For each miRNA, 4 ng of total RNA was used as 
input and miRNAs were converted to cDNA using the 
TaqMan miRNA Reverse Transcription cDNA Synthesis 
kit and miRNA-specific looped primers from the Early 
Access miRNA Profiling Kit (both from Applied Biosys- 
tems, Foster City, CA) following the manufacturer's 
recommendations. Real-time quantitative PCR was per- 
formed using miRNA-specific primer/probe pairs from the 
Early Access miRNA Profiling Kit and reactions were car- 
ried out as described elsewhere [36]. Actual amplification 
and data collection were performed on the ABI 7700 
Sequence Detection System (Applied Biosystems, Foster 
City, CA). After visual inspection, data were exported to a 
text file and further analyzed in Microsoft Excel (Microsoft 
Corporation, Redmond, WA) and various Bioconductor 
packages [37] in the R statistics environment [38]. 

Statistics for miRNAs expression data 

Delta Ct values for each miRNA were calculated with 
/2sa-let-7a as a reference, according to the manufacturer's 
directions. Only miRNAs with a dCt < 12 (where dCt = 
Cttarget " Ct/,sa_iet-7a) Were Considered for further analysis. 
Data were combined into a convenient ExpressionSet 
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(eSet) structure using the Biobase package, and further 
analyzed by means of gene-by-gene one-factor ANOVA 
using the package LMGene, all in Bioconductor and R. 
Only genes that had a False Discovery Rate-adjusted 
p-value of < 0.05 were considered differentially expressed. 
These genes were further assessed using non-parametric 
pairwise comparison of the different myeloid cell subsets 
using Tukey's post-hoc test in the R statistics package. 
Genes with a p-value < 0.05 were considered differen- 
tially expressed. 

Determination of TFBSs in microRNA promoters 

Putative promoter regions extending 2 kb upstream of 
the miRNAs were extracted from the Genome Browser 
sno/miRNA track of the UCSC March 2006 human gen- 
ome assembly [39,40]. The program Clover [41] was used 
to screen for over-represented TFBSs in these sequences 
using a precompiled library of TFBS motifs. The library 
contained 263 TFBS motifs (position-specific weight 
matrices) constructed from the JASPAR core database 
(2005) [42] and TRANSFAC version 7.0 [43]. 

To determine over-represented TFBS motifs. Clover 
starts for every location in a DNA sequence, by calculating 
a score reflecting the likelihood that a certain TF binds at 
that location. This score is a likelihood ratio, with in the 
numerator the probability that the sequence matches the 
positional weight matrix of the motif (the product of the 
frequencies of the nucleotides in the weight matrix at the 
positions corresponding to those in the sequence) and in 
the denominator the likelihood that the sequence is 
derived from a random sequence (the product of the fre- 
quencies of the nucleotides in each position in the back- 
ground sequence). This likelihood ratio is then averaged 
per complete sequence and over all subsets of the set of 
sequences. Finally a "raw score" is derived, by taking the 
logarithm of the averaged likelihood ratio. A raw score 
above zero thus signifies over-represented motifs and a 
raw score below zero signifies under-represented motifs. 
The raw score increases when more of the sequences con- 
tain good motif matches, and also when there are more 
good matches per sequence. The p-values for over-repre- 
sented motifs are subsequently derived by analyzing, 
whether random sets of promoter sequences from all the 
genes in the genome are likely to have the same, or a 
higher "raw score". Clover examines the over-representa- 
tion of multiple TFBSs and takes into account the multiple 
testing issue for its p-value calculation. Besides the raw 
score, Clover also reports an instance score, which is the 
logarithm of the likelihood ratio mentioned above for any 
specific TFBS at any specific position in the sequence. 

Our inputs to Clover are the miRNA promoter 
sequences and the 263 TFBS motifs mentioned above. 
For statistical calculations, promoter regions 2 kb 
upstream of all human genes in the genome are included 



as background. Our thresholds for Clover outputs are 
instance score > 6 for recognizing a specific TFBS, and 
p-value < 0.05 for over-represented motifs (default values 
in Clover). In addition, any mention of a TF motif score in 
subsequent sections refers to the instance score for the 
TFBS at a specific location in a promoter sequence. Our 
motivation for including instance scores is to enable us to 
calculate the conservation of the nucleotides at high-scor- 
ing TFBSs. 

Calculation of evolutionary conservation score of TFBSs 

In order to take into account conservation of TFBSs in the 
miRNA promoters, the PhastCons conservation track 
from the UCSC Genome Browser of January 2009 (http:// 
hgdownload.cse.ucsc.edu/goldenPath/hgl8/phastCons44- 
way/) was used, which represents the conservation of each 
nucleotide across 44 placental mammals as calculated 
using the PhastCons program [33]. The mammalian con- 
servation track was used because we examined the expres- 
sion of mammalian miRNAs that are generally not 
conserved outside this clade [44,45]. The base-by-base 
conservation scores are derived from a two-state phylo- 
HMM and are defined as the posterior probability that the 
corresponding alignment column was generated by the 
conserved state and not the non-conserved state of the 
phylo-HMM used in the calculation [33]. The score ranges 
from 0 to 1; the higher the value, the more conserved the 
nucleotide. Using in-house Python scripts, we obtained 
the conservation score for a TFBS as the median of the 
PhastCons scores of the bases in that TFBS. 

Microarray sample preparation of DCs and microarray 
analysis 

From one donor, 3 different samples (1 x 10^ cells; tech- 
nical replicate) were generated per DC subtype (imma- 
ture monocyte-derived DCs, maturing DCs, tolerogenic 
DCs and activated tolerogenic DCs) and RNA for micro- 
array analysis was isolated using the RNeasy Total RNA 
Extraction kit (Qiagen, Venlo, the Netherlands). Quality 
control, conversion to labeled RNA, hybridization and 
scanning were all performed at the microarray facility of 
the Department of Human Genetics (Radboud University 
Nijmegen Medical Centre Nijmegen, Nijmegen, the 
Netherlands), using Affymetrix technology while follow- 
ing the manufacturer's protocols. The .CEL files were 
processed using the Bioconductor package limma in the 
R statistics environment and further analyzed with the 
package panp to generate presence/absence calls from 
the microarray data. In this study, intensities above the 
p-value cut-off of 0.01 indicated presence, p-values 
between 0.01 and 0.02 indicated marginal presence, and 
p-values above 0.02 indicated absence. Complete, 
MIAME-compliant datasets were deposited with the 
Gene Expression Omnibus of the National Center for 
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Biotechnology Information (http://www.ncbi.nlm.nih. 
gov/geo/) and can be accessed through GEO Series 
accession number GSE23371. 

Analysis of miRNA targets 

In addition to identifying miRNAs that are expressed in 
DCs and monocytes, it is important to identify some of 
their target genes. This may help in understanding the 
miRNA-related molecular mechanisms underpinning DC 
maturation. We used TargetScanHuman version 5.1 [15] 
to predict the target genes for the miRNAs identified to 
be expressed in DCs and monocytes. The precompiled 
predictions were downloaded from the official TargetS- 
can website (http://www.targetscan.org/) and were 
further filtered for RefSeq transcripts having at least one 
site that is conserved across placental mammals and also 
had a target score (total context score) of < -0.4 before 
further analyses. To improve the reliability of the target 
genes predicted, we used a second miRNA target predic- 
tion tool, PicTar [46]. The precompiled miRNA target 
predictions were downloaded from the PicTar website at 
http://pictar.mdc-berlin.de, on January 2011. We filtered 
the predictions at a PicTar threshold of >0.4. For further 
analysis we used the target genes that were predicted by 
both TargetScan and PicTar and also by each of the algo- 
rithms separately. To narrow the scope of the target 
genes to DCs, we used the previously published dataset 
of Lehtonen and colleagues wherein the relative expres- 
sion levels of gene expression in DCs relative to mono- 
cytes are reported [47]. Using their data, we selected 
genes that are at least 2-fold down-regulated in DCs rela- 
tive to monocytes and that we also predicted to be target 
genes for DC-expressed miRNAs. Comparisons between 
our predicted target genes and genes that were over- or 
under-expressed in DCs relative to monocytes in the 
Lehtonen dataset were done at the level of RefSeq DNA 
ID to ensure that the correct gene isoforms were being 
matched. (The Lehtonen dataset was based on global 
gene expression analyses using Affymetrix HG-U133A 
Gene Chip oligonucleotide arrays. We used Biomart in 
Ensembl version 62 available at http://www.ensembl.org/ 
biomart to obtain the corresponding RefSeq DNA IDs 
for the probesets in the microarray dataset). The set of 
genes so selected are likely to be important in miRNA 
regulation of DC maturation from monocytes. 

Gene Ontology analysis 

Gene ontology enrichment analyses were done using the 
Cytoscape [48] plugin BINGO version 2.3 [49]. Using this 
software, we tested for over-representation using the 
hypergeometric test with Benjamini & Hochberg False 
Discovery Rate (FDR) correction. GO processes reported 
were deemed significant when the corrected p-value was < 
0.05. Venn diagrams were drawn using the Venny tool of 



Oliveros, J.C. available at http://bioinfogp.cnb.csic.es/tools/ 
venny/ index.html. 

Results 

Monocytes and monocyte-derived dendritic cell subtypes 
express distinct microRNAs 

Monocytes were isolated to > 90% purity and monocyte- 
derived dendritic cell subtypes were generated and 
validated as depicted in Additional file 1, Figure SI. A 
TaqMan-based quantitative RT-PCR method was used to 
profile the expression of miRNAs. Immature monocyte- 
derived DCs (iDCs) were induced to either undergo 
maturation by LPS stimulation for 6 hr (mDCs), to mature 
into tolerogenic DCs (tDCs) in the presence of IL-10 and 
dexamethasone for 24 hr, or to mature into activated tol- 
erogenic DCs (atDCs; previously described by Emmer et 
al. [50]) in the presence of IL-10 and dexamethasone for 
24 hr, followed by 6 hr of LPS (Additional file 1, Figure 
SIA). Phenotypical and functional characterization con- 
firmed that mature DCs expressed high levels of CD80 
and CD86 (Additional file 1, Figure SIB), induced strong 
proliferation of allogeneic T cells in a mixed lymphocyte 
reaction (Additional file 1, Figure SIC) and secreted high 
levels of pro-inflammatory cytokines TNF and IL8 (Addi- 
tional file 1, Figure SID). Accordingly, iDCs and tDCs did 
not show up-regulation of CD80 and CD86, induced 
much less proliferation and produced less cytokines. Acti- 
vated tolerogenic DCs showed a more intermediate pat- 
tern of activation, indicative of less profound TLR 
signaling in response to LPS. The expression profile of 157 
miRNAs in monocytes and the various DC activation 
stages was determined in a primary screen involving 
total RNA from one donor as shown in Additional file 2, 
Figure S2. 

Of the 157 miRNAs screened, 104 were detected and 
expressed to variable degrees in monocytes and/or DCs, 
whereas 53 miRNAs remained undetected, irrespective 
of cell type or stimulus (Additional file 3, Table SI). 
Intriguingly, 27 miRNAs appeared to be differentially 
expressed between monocytes, DCs and stimulated DCs 
(Table 1). These 27 miRNAs were selected for additional 
profiling of multiple donors as described below. 

Identification of miRNAs differentially regulated during 
monocyte-derived DC activation 

To validate the expression of the identified 27 miRNAs, 
their expression was further examined in monocytes and 
DCs of 3 different donors (Additional file 2, Figure S2). 
Statistical analysis revealed that 18 miRNAs were truly 
differentially expressed (gene-by-gene ANOVA, p < 0.05; 
Table 2). Six miRNAs were up-regulated in monocytes 
versus DCs, and 10 were up-regulated in DCs versus 
monocytes. Eight miRNAs were not differentially 
expressed and had a high degree of variation in 
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Table 1 MicroRNAs rescreened after primary screen in monocytes and dendritic cells, and their differences in 
expression level between cell types 



MicroRNA in 


Mature miRNA Sequence 


Current genomic entries in miRbase 


Up in 


Up in DC Up in No 


assay 






matching mature miRNA 


DC 


subset monocytes difference 


hsa-let-7e 


UGAGGUAGGAGGUUGUAUAGU 


same 


X 




teo-mlR-' 


iSb 


UAGCAGCACAUCAUGGUUUACA 


same 




X 


/isa-mlR-' 


16 


UAGCAGCACGUAAAUAUUGGCG 


hsa-miR-1 6-1 , hsa-miR-16-2 




X 


hsa-miR-27a 


UUCACAGUGGCUAAGUUCCGCC 


same 




X 


/i5a-miR-27b 


UUCACAGUGGCUAAGUUCUG 


same 




X 


feo-miR-34a 


UGGCAGUGUCUUAGCUGGUUGU 


same 


X 




/isa-miR-99a 


AACCCGUAGAUCCGAUCUUGUG 


same 


X 




hsa-m\R-' 


100 


AACCCGUAGAUCCGAACUUGUG 


same 


X 




/isa-mlR-' 


125a 


UCCCUGAGACCGUUUAACCUGUG 


same 


X 




hsa-mlR-' 


125b 


UCCCUGAGACCGUAACUUGUGA 


hsa-miR-1 25b-l, hsa-mlR-125b-2 


X 




/isa-mlR-' 


126 


UCGUACCGUGAGUAAUAAUGC 


same 




X 


hsa-m\R-' 


130a 


CAGUGCAAUGUUAAAAGGGC 


same 




X 


/iso-mlR-' 


132 


UAACAGUCUACAGCCAUGGUCG 


same 




X 


hsa-mlR-' 


135a 


UAUGGCUUUUUAUUCCUAUGUGA 


hsa-miR-135a-1, hsa-mlR-135a-2 


X 




hsa-m\R-' 


137 


UAUUGCUUAAGAAUACGCGUAG 


same 


X 




hsa-m\R-' 


140 


AGUGGUUUUACCCUAUGGUAG 


same 




X 


/isa-mlR-' 


146 


UGAGAAGUGAAUUCCAUGGGUU 


hsa-miR-146a 




X 


hsa-m\R-' 


150 


UCUCCCAACCGUUGUACCAGUG 


same 




X 


/isa-mlR-' 


155 


U U AAUGCUAAUCGUGAUAGGGG 


same 




X 


feo-mlR-' 


186 


CAAAGAAU UCUCCU U U UGGGCU U 


same 




X 


/isa-mlR-' 


199b 


CCCAGUGUUUAGACUAUCUGUUC 


same 




X 


/isa-mlR-' 


199-s 


CCCAGUGUUCAGACUACCUGUU 


hsa-miR-1 99a-1, hsa-mlR-199a-2 




X 


fca-miR-210 


CUGUGCGUGUGACAGCGGCUG 


same 




X 



expression levels (data not shown). However, hsa-miR- 
15b and /jsa-miR-16 showed remarkably similar expres- 
sion profiles and are expressed at higher levels in mono- 
cytes (Figure 1, upper row of charts; Pearson correlation 
= 0.99). Hsa-miR-125a, fafl-miR-221 and fofl-miR-342 
are examples of miRNAs that are up-regulated in the DC 
populations relative to monocytes (Figure 1, middle row). 
Furthermore, both /2sa-miR-146a and fafl-miR-155 were 
up-regulated during differentiation, and were slightly up- 
regulated in mDCs and atDCs (Figure 1, two left charts 
in lower row; gene-by-gene ANOVA, p < 0.05, Tukey's 
post-hoc test; Pearson correlation = 0.90), whereas hsa- 
miR-210 was primarily up-regulated in atDCs (Figure 1, 
right chart of lower row). Data generated from these ana- 
lyses indicate that many of the selected miRNAs are 
indeed expressed in monocytes and DCs, and that 3 of 
these miRNAs are differentially regulated during DC 
development and maturation. 

Identification of TFBSs In promoters of differentially 
expressed miRNAs in monocytes and DCs 

In order to link the miRNA expression profiles to a spe- 
cific TFBS (or groups of TFBSs) of the myeloid cell 



types that were analyzed, sequences 2 kb upstream of 
the mapped pre-miRNAs were chosen for analysis. It 
should be noted that it is expected that a majority of 
miRNAs contain their own promoters [21], and that in 
some cases, different miRNA genes (such as hsa-miR- 
16-1 and /zsfl-miR-16-2) in the genome give rise to an 
identical mature miRNA (such as hsa-miR-16) . There- 
fore 31 promoters could be assigned to 27 miRNAs that 
were profiled in detail (Table 1). Likewise, 167 promo- 
ters could be assigned to the total of 157 miRNAs ana- 
lyzed. Analysis of the 31 promoters revealed that 3 
TFBS motifs (Elk-1, RREB-1 and SPIB) were over-repre- 
sented using the Clover criteria (p < 0.05) among all 8 
promoter sequences of miRNA genes that are up-regu- 
lated in monocytes, and to have at least 1 high-scoring 
TFBS (instance score > 6) in each promoter. Selecting 
TFBSs using these same criteria (over-represented 
motifs and at least 1 high-scoring TFBS per promoter) 
among the 12 promoters of the miRNAs that are up- 
regulated in DCs relative to monocytes, reveals that they 
have 13 TFBSs in common (Figure 2, Table 3, Addi- 
tional file 4, Table S2). Furthermore, 4 TFBSs are over- 
represented in the promoters in all of the 8 miRNAs 
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Table 2 Assignment of p-values for expressed miRNAs in 
all cell types (monocytes and various DC subsets) and DC 
subsets alone after ANOVA 



Posterior p-values* 


miRNA gene 


All cell types 


DC subsets 


hsa-let-7e 


1 .9E-05 


4.2E-01 


hsa-miR-15b 


1 .2E-04 


6.8E-01 


hsa-miR-16 


2.1E-02 


7.8E-01 


hsa-miR-27a 


3.0E-01 


2.5E-01 


hsa-miR-27b 


8.2E-01 


8.2E-01 


hsa-miR-34a 


3.0E-05 


9.6E-01 


hsa-miR-99a 


2.2E-06 


7.2E-01 


hsa-miR-100 


3.7E-06 


73E-01 


hsa-miR-125a 


1 .3E-04 


8.6E-01 


hsa-miR-125b 


1 .6E-05 


7.2E-01 


hsa-miR-126 


6.1E-01 


7.7E-01 


hsa-miR-130a 


1 .6E-03 


44E-01 


•*hsa-miR-132 


4.0E-02 


9.8E-01 


hsa-miR-135a 


7.1E-10 


9.1E-02 


hsa-miR-137 


3.2E-07 


2.1E-01 


hsa-miR-140 


8.2E-01 


7.1E-01 


1 lia 1 1 1 1 r\ 1 T^ua 


9 6E-05 


1 3E-02 


•*»hsa-miR-150 


9.9E-02 


7.1E-01 


hsa-miR-155 


4.7E-06 


2.5E-02 


hsa-miR-186 


7.2E-01 


6.9E-01 


hsa-miR-199b 


5.2E-05 


7.1E-01 


hsa-miR-199s 


84E-03 


44E-01 


hsa-miR-210 


1 .3E-05 


4.0E-02 


hsa-miR-221 


5.6E-06 


74E-01 


hsa-miR-326 


74E-01 


9.7E-01 


hsa-miR-340 


8.6E-01 


7.9E-01 


hsa-miR-342 


6.1E-06 


9.9E-01 



*) Posterior p-values < 0.05 indicate differential expression among subsets. 
**) Considered not expressed differentially after visual inspection of gene 
expression data. 

***) Considered expressed differentially after visual inspection of the data. 



that did not show any appreciable difference in expres- 
sion levels (Table 3). The complete set of high-scoring 
TFBSs found in the promoters of the miRNAs that are 
up-regulated in DCs and monocytes are summarized in 
Figure 3 and Additional file 5, Figure S3 respectively. 
Clover detects motifs that are over-represented relative 
to their expected frequencies and calculates p-values 
based on a comparison with all promoters in the gen- 
ome. Nevertheless such predictions can potentially give 
rise to false positives, as even random sets of miRNA 
promoters contain some over-represented motifs (Figure 
3B). We therefore also examined the over-representation 
of motifs in the promoters of DC and monocyte miR- 
NAs relative to those of random sets of miRNAs, 



selected from the total set of 167 miRNA promoters. 
The 12 miRNAs that are over-expressed in DCs share 
more TFBS motifs than do random sets of 12 miRNA 
promoters, specifically with respect to motifs that are 
shared by many of the promoters (Figure 3B). When 
quantifying the number of shared motifs as the sum of 
all occurrences of all Clover detected motifs, across all 
promoters, we observe an enrichment of 1.96 relative to 
randomly selected miRNA promoters (p = 0.001). Like- 
wise, using random sets of 8 miRNA promoters, we 
observed that the 8 miRNAs that are over-expressed in 
monocytes have a factor of 1.58 enrichment of TFBS 
motifs (p = 0.046; Additional file 5, Figure S3B). The p- 
value here was estimated as the number of times, out of 
1000, when a randomly chosen set of miRNA promoters 
had at least an equal number of motifs as did the test 
set of sequences. Given the significant over-representa- 
tion of different sets of TFBSs in the different sets of 
miRNAs, it appears that miRNAs contain multiple, dif- 
ferent TFBSs in their promoters that specify their 
expression in certain cell types. This is in agreement 
with other studies that show that transcription factors 
often work in combinations [25,26]. 

Evolutionary conservation and distribution of TFBSs in 
the promoter sequences of myeloid miRNAs 

Assuming that regulation of miRNA expression is con- 
served among mammals, we used the PhastCons evolu- 
tionary conservation track that is based on nucleotide 
conservation among placental mammals (http://hgdown- 
load.cse.ucsc.edu/goldenPath/hgl8/phastCons44way/) to 
quantify evolutionary conservation of the predicted 
TFBSs. We compared both the motif instance scores 
and raw scores of TFBSs, as derived from Clover, with 
the nucleotide conservation score of the TFBSs. Applied 
to the 2 kb upstream sequences of the pre-miRNAs that 
are up-regulated in DCs, a small but significant correla- 
tion of 0.12 (p = 6.17 X 10'^) was observed between the 
conservation scores and the TFBS instance scores. Simi- 
lar results were obtained using the raw score (r = 0.22, 
p = 8.99 X 10'^). An overlay of TFBS instance scores 
and conservation scores for the DC miRNA promoters 
is depicted in Figure 4A. The positive but weak correla- 
tion between TFBS motif scores and conservation scores 
suggests that they are largely independent measures and 
as such, combination of the two as a filtration measure 
would be non-redundant. When filtering predicted sites 
based both on their conservation and on their TFBS 
instance scores, TFBSs tend to cluster in a region 500 
bp just upstream of the putative TSS of miRNA genes 
that are over-expressed in DCs (Figure 4B). This result 
corroborates other findings indicating that TFBSs tend 
to cluster near the transcription start site (TSS) of genes 
[51-54]. In addition, it suggests that combining the 
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motif matching score with the extent of evolutionary 
conservation of a TFBS would likely reveal TFs biologi- 
cally relevant to the identified set of miRNA promoters. 

In further support of the conservation of the binding 
motifs that are over-represented in the promoters of the 
DC-expressed miRNAs, we created a set of random 
sequence blocks of the same length-distribution as the 
predicted TFBSs from the same promoters as the test 
set. From these two data sets, we observed that the pre- 
dicted TFBSs are significantly (p = 1.49 x 10 ^ Wil- 
coxon rank sum test) more conserved than randomly 
chosen sequences. The conservation score assigned to a 
TFBS is the median of the PhastCons scores of the 
bases in that TFBS segment. We observed that the med- 
ian of the conservation scores of the test TFBSs is more 
than 2 fold higher than that of the random sequences 
(Figure 4C). The complete set of conserved high-scoring 
TFBSs found in the promoters of the miRNAs up-regu- 
lated in DCs and monocytes are summarized in Figure 5 
and Additional file 6, Figure S4 respectively. Further- 
more, we observed at the threshold of 10''^ percentile of 
median TFBS conservation, that the miRNAs that are 



over-expressed in DCs share more (fold = 3.51, p < 
0.001) sites for TF motifs than did 1000 random sets of 
promoters of equal size, length and at the same conser- 
vation threshold (Figure 5B). A similar result (fold = 
3.48, p = 0.003) was obtained for the miRNAs that were 
over-expressed in monocytes (Additional file 6, Figure 
S4B). These results show that evolutionary conserved 
TFBSs are more common in the promoters of miRNAs 
over-expressed in DCs and monocytes than in randomly 
chosen miRNA promoters. 

Integrating the motif-matching and evolutionary 
conservation scores of TFBSs and further validation of 
cognate TFs in microarray data from DCs 

Without taking evolutionary information into account, 13 
TFBSs were found to occur at least once in all the promo- 
ters of miRNAs preferentially expressed in DCs (Figure 2). 
Filtering TFBSs at 10* percentile threshold of the TFBS 
conservation scores not surprisingly identifies TFBSs that 
are common to far less miRNAs promoters (Figure 5). 
Among these, the TFBS of MZFl was the most conserved 
at the thresholds used. To validate the expression of the 
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cognate transcription factors whose binding site was iden- 
tified, a transcriptome was generated of DCs from one 
donor using microarrays. The resulting data were screened 
for the presence or absence of transcription factors (TFs) 
whose binding sites were predicted in promoters of miR- 
NAs up-regulated in DCs. For 10 out of the 17 TFBSs that 
are over-represented among the DC promoters (Addi- 
tional file 4, Table S2B) and that contained at least one 
TFBS with an instance score of at least 6.0 and a conserva- 
tion threshold score of 0.968 (i.e. the 10* percentile of the 
median conservation scores of all TFBSs predicted in the 
set of miRNA promoters) as shown in Figure 5, the 
mRNA encoding their cognate TF was expressed in DCs. 
Expression of this 58.8% of the predicted over-represented 
TFs is slightly higher than that for all the TFs that were 
examined with the microarrays (54% expressed, 46% not 
expressed. Additional file 7, Table S3). Nevertheless, we 



found a high correlation (r = 0.91, p = 1.17 x 10'^) when 
the number of motifs of expressed TFs was correlated 
with the number of miRNA promoters in which the motifs 
of these TFs were detected (Additional file 8, Figure S5A). 
An insignificant correlation (r = -0.72, p = 4.88 x 10' ) 
was observed when a random set of miRNA promoters 
was used (Additional file 8, Figure S5B). 

Gene Ontology of the expressed TFs highlighted 
immune-related processes like "interspecies interaction 
between organisms" (p = 4.49 x 10'^) and "regulation of 
T-helper 2 type immune response" (p = 1.19 x 10'^, 
Additional file 9, Table S4A). This is in line with the 
expected functions of DCs in immune response. Mean- 
while, Gene Ontology of the TFs that were not 
expressed highlighted development-related processes 
like "cell development" (p = 1.88 x 10' ) and "cell fate 
commitment" (p = 7.54 x 10'^, Additional file 9, Table 



Table 3 TFBSs that are present in all miRNA promoters In a particular (group of) cell type(s) 



Cell types 


Overrepresented TFBS 


Monocytes 


Elk-1, RREB-1, SPIB 


DCs* 


AP-4, E47, MAZR, Myf MZFl, RREBl, RREB-1, Spl, SPIl, SPIB, STAT6, ZNF42_l-4, ZNF42_5-13 


Subset of DCs 


AP-1, ARP-1, FOXn, myogenin/NF-1, Ncx, Oct-1, OCT-x, Spl, TCFl 1, ZNF42_l-4 


All cells"* 


Oct-1, RP58, SPIB, USE 


*) Includes monocyte-derived iDCs, maturing DCs, tolerogenic DCs and activated tolerogenic DCs 
**) Includes maturing DCs and/or activated tolerogenic DCs 

***) Concerns all cells investigated, and miRNAs did not show appreciable differences in expression level 
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S4B), with little relevance to immune-related processes. 
Taken together, these data provide evidence that many 
expressed miRNAs in DCs have evolutionarily conserved 
TFBSs that may be relevant to DC biology in their 
promoters. 

Target prediction and Gene Ontology analysis of 
microRNAs enriched in monocytes and DCs 

To extend the regulatory pathways in which the miRNAs 
participate downstream, we predicted miRNA targets, 
using TargetScanHuman version 5.1 [15] and PicTar 



[46]. Both datasets were based on conservation in mam- 
mals. We obtained the precompiled dataset from their 
official websites. At a TargetScan context score threshold 
of at most -0.4, we obtained 421 distinct target genes 
(Additional file 10, Table S5; counting in RefSeq DNA 
IDs) for the miRNAs that are over-expressed in DCs. 
Meanwhile at a PicTar score threshold of at least 0.4, we 
obtained 2300 distinct target genes (Additional file 11, 
Table S6) for the miRNAs over-expressed in DCs. Like- 
wise, we predicted 282 (Additional file 12, Table S7) and 
1591 (Additional file 13, Table S8) distinct genes using 
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TargetScan and PicTar respectively, to be targets for the 
miRNAs that are over-expressed in monocytes. For both 
DC-over-expressed miRNAs and monocyte-over- 
expressed miRNAs, the overlap of targets predicted by 
both TargetScan and PicTar is less than 50% of their 
respective total predictions (Additional file 14, Figure S6) 
suggesting that we still need to examine their indepen- 
dent predictions. Interestingly, only 7 identical targets 
were predicted for both monocyte-over-expressed 



miRNAs and DC-over-expressed miRNAs at the TargetS- 
can threshold of <-0.4 and PicTar threshold of >0.4. The 
numbers of target genes predicted are summarized in 
Additional file 14, Figure S6. 

The number of target genes predicted using TargetS- 
can and PicTar are very large even at the stringent 
thresholds used. Furthermore, Gene Ontology (GO) ana- 
lyses of the gene sets did not yield conclusive results. To 
increase the relevance of the predicted targets to DCs 
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we further filtered them based on their differential 
expression in DCs relative to monocytes during DC dif- 
ferentiation using previously published data by Lehtonen 
et al (2007) [47]. At a gene expression fold change 
threshold of 2, we identified 34 unique genes corre- 
sponding to the DC-associated targets predicted using 
PicTar to be down-regulated at all time points in DCs 
relative to monocytes (Additional file 15, Table S9A). 
Amongst this filtered set of genes is ILIO, which inter- 
estingly is involved in the process of "negative regulation 
of myeloid dendritic cell activation" (p = 2.64 x 10'^, 
from GO analysis of the 34 down-regulated genes). Cor- 
respondingly, 10 unique genes were associated with the 
down-regulated genes of the TargetScan-predicted tar- 
gets of miRNAs that were over-expressed in DCs rela- 
tive to monocytes (Additional file 15, Table S9B). In 
total, 6 of these down-regulated genes [SlKl, PELI2, 
MYLIP, RNF138, CCNG2 and FOXOl; in order of 
increasing down-regulation of their RNA transcripts at 
time point 24 hr) were also predicted using PicTar 
(Additional file 15, Table S9). Of these, SIKl, PELI2 and 
one transcript of CCNG2 were increasingly down-regu- 
lated in DCs relative to monocytes in the time course of 
DC differentiation from monocytes; suggesting their 
involvement in mechanisms of DC differentiation. 

In addition to the DC dataset, we also examined the 
expression of the predicted target genes of the miRNAs 
that are over-expressed in monocytes relative to DCs. We 
found 35 unique genes associated with the PicTar-pre- 
dicted targets of miRNAs over-expressed in monocytes 
relative to DCs that were down-regulated in monocytes 
relative to DCs (Additional file 16, Table SlOA). Corre- 
spondingly, transcripts of 9 unique genes associated with 
the TargetScan-predicted targets of miRNA over- 
expressed in monocytes relative to DCs were also down- 
regulated at all time points in monocytes relative to DCs 
(Table SlOB). Among these is DICERl which is involved 
in RNA interference and is up-regulated in DCs. Further- 
more, there was an overlap of 2 genes {ACVR2A and 
MAP3K4) between the PicTar and TargetScan target gene 
sets (Additional file 14, Figure S6 and Additional file 17, 
Table Sll). One transcript o(ACVR2A was increasingly 
down-regulated (Additional file 16, Table SIO). Moreover, 
there was no overlap in expressed genes that were targets 
to both DC- and monocyte-over-expressed miRNAs, sug- 
gesting congruency in our data. These results suggest that 
miRNA target genes can be either up-regulated or down- 
regulated in myeloid cells to regulate differentiation of the 
myeloid cells. 

Discussion 

As little is known about the expression and regulation 
of miRNAs in monocytes, DCs and their stimulated pro- 
geny, part of the miRNA transcriptome of these cells 



was generated and analyzed. Out of 157 miRNAs pro- 
filed, 104 appeared to be expressed in myeloid cells to a 
varying degree. Since the database miRBase 13.0 (March 
2009) reveals the existence of at least 706 human miR- 
NAs, it is conceivable that a much higher number of 
miRNAs is expressed in human monocytes and DCs. Sbc 
miRNAs appear up-regulated in monocytes, 10 are up- 
regulated in DCs and 3 are differentially expressed 
among different DC populations in response to LPS. 
Also, some miRNAs showed similar changes in expres- 
sion levels among DC subsets. 

Amongst the 10 miRNAs that were over-expressed in 
DCs as compared to the monocytes from which they 
were derived, /z5fl-miR-34a, fafl-miR-125a, fafl-miR-342 
and hsa-let-7e have been shown to be up-regulated in 
the course of DC differentiation from monocytes in cul- 
ture [55,56]. The up-regulation of /zsfl-miR-342 in 
monocyte-derived DCs as compared to monocytes is 
likely the result of culture conditions, as fafl-miR-342 is 
expressed at a much lower level in freshly isolated 
blood-derived myeloid DCs than in DCs generated in 
vitro (data not shown). Of the remaining miRNAs, 6 
show an upward trend in expression level in blood- 
derived DCs, but the levels in cultured DCs are higher. 
Intriguingly, /jsa-miR-146a and hsa-miR-155 appear to 
respond to LPS, but much less so than in monocytes 
and macrophages [7,8]. Recently, it was demonstrated 
that LPS-mediated activation of protein kinase Aktl 
results in up-regulation of miRNA let-7e in primary 
macrophages, while at the same time repressing miR- 
155 expression [57]. Our data show that both hsa-let-7e 
and hsa-miR-155 are up-regulated in DCs compared to 
monocytes, but that only hsa-miR-155 is slightly up- 
regulated by LPS. Taken together, these data imply that 
intrinsic differences between DCs and macrophages 
exist in response to the TLR4 ligand LPS. 

With regard to the targets of the miRNAs screened 
herein, recent literature indicates that inhibition of hsa- 
miR-34a or addition of one of its target genes, JAGl, 
have been observed to functionally stall the differentia- 
tion of monocyte-derived dendritic cells [55]. This sup- 
ports our identification of /zsfl-miR-34a as an over- 
expressed miRNA in DCs relative to monocytes and 
JAGl as its target gene (NM_000214 in Additional file 
11, Table S6). In addition, we provide a list of predicted 
DC targets that are down-regulated in expression in 
DCs relative to monocytes (Additional file 15, Table S9). 
In this list, one of the targets of the DC-expressed 
miRNA, fafl-let-7e, is ILIO, which is involved in negative 
regulation of myeloid dendritic cell activation. Taken 
together, these data highlight miRNAs, and their target 
genes, that can potentially modulate DC differentiation. 

Extensive miRNA promoter analysis revealed that 13 
TFBSs are over-represented and commonly shared by 
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the 12 promoter sequences of the 10 miRNAs that are 
up-regulated in DCs. These include signal transducer 
and activator of transcription 6 (STAT6), which is 
known to be involved in IL4 signaling and DC differen- 
tiation and maturation [58], spleen focus forming virus 
(SFFV) proviral integration oncogene spil (SPIl, also 
known as PU.l), which is indispensable for normal mye- 
loid and lymphoid development [59] and specificity pro- 
tein 1 (SPl), which is involved in the expression of 
dendritic cell-specific ICAM-3 grabbing non-integrin 
(DC-SIGN) [60]. 

Although the TSSs in these promoter regions are not 
known, TFBSs did cluster, especially when taking evolu- 
tionary conservation into account, within the 500 bp 
upstream region of the annotated pre-miRNAs, coincid- 
ing with the region in miRNA promoters at which TSS 
have been discovered experimentally [19,61]. Amongst 
the TFBSs predicted in the DC miRNAs, the sites for 
myeloid zinc finger (MZFl) is best conserved of all, even 
though the role of its cognate TF in DCs remains elu- 
sive. Nevertheless, MZFl is thought to be a bi-functional 
transcriptional regulator, repressing transcription in 
non-hematopoietic cells, activating transcription in cells 
of hematopoietic origins and controlling cell prolifera- 
tion and tumorigenesis [62,63]. In addition, one of the 
conserved TFBSs for which the cognate TF was not 
expressed in monocyte-derived DCs appears to be Spi-B 
transcription factor (SPIB), which has been implicated in 
plasmacytoid DC (pDC) development, a non-myeloid 
cell type [64,65]. It should be noted, however, that many 
of the miRNAs expressed in myeloid monocyte-derived 
DCs are also expressed in pDCs (data not shown). 

When taking all predicted and well-conserved TFBSs 
in the promoters of miRNAs that are up-regulated in 
DCs into account, the number of TFs expressed 
increases with the degree at which their TFBS motifs 
are shared between the promoter sequences. Gene 
Ontology analysis indicates that the type of expressed 
TFs enriched in these data are relevant to the immune 
system process and are in line with the known function 
of DCs, whereas those that were not expressed are rele- 
vant to cell development (Additional file 9, Table S4). 
Importantly, libraries of TFBS motifs used do not repre- 
sent all possible TFBSs and as such not all possible 
TFBSs for expressed miRNAs have been identified in 
this study. Moreover, for consistency, we have used 
TFBS motifs, and not TFs, in making comparisons of 
TFBSs because the databases used have redundant motif 
names for the same TFs. It should be noted that pre- 
dicted target genes were not uniformly down-regulated, 
as we found evidence in DCs whereby predicted target 
genes were up-regulated in the cells or were also targets 
for the miRNAs that were over-expressed in monocytes. 
This may be due to limitations of the target prediction 



algorithms, or the targets might still be down-regulated 
at the protein level. Nevertheless, none of the down- 
regulated DC-miRNA targets was a predicted target for 
the monocyte-miRNAs. Furthermore, there is evidence 
that in a minority of cases, target genes are actually up- 
regulated by miRNAs [66,67]. The up-regulation of 
DICERl in DCs relative to monocytes is of special inter- 
est, as it is known that innate immune signaling is 
tightly controlled by miRNAs [7,8]. Furthermore, there 
is also evidence that one type of DC, the Langerhans 
cell, requires proper functioning of DICERl to induce 
CD4 T cell function [68]. Together with the fact that 
there are more up-regulated miRNAs in DCs than 
monocytes in our data set, it is tempting to speculate 
that DICERl up-regulation is required for proper DC 
function. 

Conclusions 

The data provide evidence that, among the many 
expressed miRNAs in DCs, evolutionarily conserved 
TFBSs relevant to DC biology are present in their pro- 
moters. Furthermore, the identified miRNAs, their asso- 
ciated TFs and predicted target genes could help 
improve our understanding of the molecular pathways 
that underpin DC differentiation and maturation. 
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Additional file 1: Figure SI. Setup of culture and quality of monocyte- 
derived DCs. (A) Culture setup and harvest schedule for RNA isolations. 
(B) Purity of monocytes (upper histogram) and the expression of 
maturation markers CD80 and CD86 on the different DC populations (C) 
Mixed lymphocyte reaction witU the various DC populations; ratio of DCs 
vs. PBLs is indicated on the x-axis. (D) Production of TNFa and IL8 by the 
different DC populations, as determined by ELISA. The white and black 
bars each represent data from two different donors. 

Additional file 2: Figure S2. Schematic overview of the miRNA 
expression screen in monocytes and DCs. 

Additional file 3: Table SI. IVlicroRNAs for which assays have been 
developed by Applied Biosystems, and of which expression levels were 
determined in monocytes and dendritic cells. Of a total of 157 miRNAs, 
104 were expressed in DC and/or monocytes, whereas 53 were not 
detected in either cell type. 

Additional file 4: Table S2. Overview of over-represented TFBS motifs in 
promoter seqeunces of miRNAs that are upregulated in monocytes or 
DCs. 

Additional file 5: Figure S3. (A) Commonality of TFBSs in promoter 
regions of monocyte-expressed miRNAs. Shown on the x-axis are the 
high-scoring TFBSs (i.e. of instance score > 6) that occur at least once in 
the 2 kb promoters of the miRNA up-regulated in monocytes. The y-axis 
shows the number of promoters that have the TFBS, at least once at this 
threshold. (B) The distribution of the number of common TFBS hits per 
number of common miRNA promoters as in A, for test and random sets 
of miRNA promoters. The values for random are the median values from 
1000 random set of miRNA promoters of same size and length as those 
in the monocyte set. The score is the ratio of the sum of "the product of 
numbers on the x-axis and corresponding y-axis values" for the 
monocyte set relative to that of the random. The p-value is the fraction 
of cases wherein this sum for random sets of miRNA promoters is 
greater than or equal to that of the DC set. 
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Additional file 6: Figure S4. (A) TFBSs shared among the promoter 
regiorns of miRNAs up-regulated in monocytes. Shown on the x-axis are 
high-scoring TFBSs filtered at a motif instance score threshold of atleast 
6 and at a 10* percentile evolutionary conservation score that occur at 
least once in the promoters of miRNAs up-regulated in DCs. The y-axis 
shows the number of promoters that have the TFBS at least once at 
these thresholds. (B) The distribution of the number of common TFBS 
hits per number of common mlRNA promoters as in A, for test and 
random sets of miRNA promoters. The values for random are the median 
values from 1000 random set of miRNA promoters of same size and 
length as those in the monocyte set. The score is the ratio of the sum of 
all occurrences of all TFBS for the monocyte set relative to that of the 
random sets. The p-value is estimated from the number of instances 
wherein this sum of "the product of number of TFBSs and the number 
of miRNA promoters" of 1000 random sets of miRNA promoters is 
greater than or equal to that of the monocyte set 

Additional file 7: Table S3. Library of TFBS motif for which the cognate 
TFs are expressed (presence = 1) or not expressed (absence = 1) in DCs. 

Additional file 8: Figure S5. Correlation of the number of expressed TFs 
with the number of promoters of miRNAs that are over-expressed in 
DCs. The TFs used are those of high-scoring TFBSs (i.e. of instance score 
> = 6) that are also highly conserved (10th percentile of PhastCons 
scores) and occur at least once in the 2 kb promoters of the miRNA up- 
regulated in DCs (left plot), and a random set of miRNAs (right plot). The 
original number of miRNA promoters that share at least one TFBS was 12 
for both DCs and the random set of promoters (i.e the number of 
promoters of miRNAs that were over-expressed in DCs). Due to the high 
conservation threshold used, the maximum number of miRNA promoters 
that share at least one TFBS became smaller (6 for the test set and 3 for 
the random set). The correlations were calculated using the Pearson 
correlation as implemented in R. 

Additional file 9: Table S4. Gene ontology analysis of the expressed 
and not-expressed TFs of predicted binding sites in the promoters of 
miRNA. 

Additional file 10: Table S5. The 421 target genes that are predicted by 
TargetScan for the DC-over-expressed miRNAs. 

Additional file 11: Table S6. The 2300 target genes that are predicted 
by PicTar for the DC-over-expressed miRNAs. 

Additional file 12: Table S7. The 282 target genes that are predicted by 
TargetScan for the IVlonocyte-over-expressed miRNAs. 

Additional file 13: Table S8. The 1591 target genes that are predicted 
by PicTar for the Monocyte-over-expressed miRNAs. 

Additional file 14: Figure S6. Venn diagram showing the intersection of 
target genes. (A) Intersection of target genes of the miRNAs that are 
over-expressed in DCs, using PicTar (DC.pic), and TargetScan (DC.tar). (B) 
Similar to A, but for monocytes. (C) Intersections of the common targets 
found by PicTar and TargetScan for DC miRNAs (DC.pic.tar), the common 
targets for monocyte miRNAs (mono.pictar), and the dataset of 
Lehtonen et al. 2007 (Leh.) of genes that are regulated in DCs relative to 
monocytes during DC differentiation. (D-F) Intersection of common 
target genes from DC miRNAs, monocyte miRNAs, and genes that were 
up-regulated (upreg) or down-regulated (downreg), in DCs relative to 
monocytes in the data set of Lehtonen et al.2007. Figure D, E, F 
represent respectively the data sets at time points 3, 6 and 24 hr of DC 
differentiation from monocytes. (Differential gene relations (DC/ 
monocyte) were selected at a fold change of at least 2. The PicTar score 
used > = 0.4, TargetScan context score < = -0.4. Comparisons were done 
at the level of RefSeq DNA ID and the elements in the intersections with 
the upreg and downreg datasets are provided in Table SI 1 with gene 
symbols attached. 

Additional file 15: Table S9. Target genes that are down-regulated in 
DCs relative to monocytes and that were predicted using PicTar or 
TargetScan, to be targets of miRNAs that are over-expressed in DCs 
relative to monocytes. 

Additional file 16: Table S10. Target genes that are down-regulated in 
monocytes relative to DCs and that were predicted using PicTar or 



TargetScan, to be targets of miRNAs that are over-expressed in 
monocytes relative to DCs. 

Additional file 17: Table S11. Elements in the Venn diagram 
intersection of Figure S6D to S6F. 
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